GENERATING STABLE MODULAR GRAPHS 



STEFANO MAGGIOLO AND NICOLA PAGANI 



Abstract. We present and prove the correctness of the program boundary, 
whose sources are available at http : //people . sissa. it/"maggiolo/boundary/ 
Given two natural numbers g and n satisfying 2g + n — 2 > 0, the program 
generates all genus g stable graphs with n unordered marked points. Each 
such graph determines the topological type of a nodal stable curve of arith- 
metic genus g with n unordered marked points. Our motivation comes from 
the fact that the boundary of the moduli space of stable genus g, n-pointed 
curves can be stratified by taking loci of curves of a fixed topological type. 



Moduli spaces of smooth algebraic curves have been defined and then compact- 
ified in algebraic geometry by Deligne and Mumford in their seminal paper ^DM]- 
A conceptually important extension of this notion in the case of pointed curves was 
introduced by Knudsen [Kj . 

The points in the boundary of the moduli spaces of pointed, nodal curves with 
finite automorphism group. These curves are called stable curves (or pointed stable 
curves). The topology of one such curve is encoded in a combinatorial object, called 
stable graph. The boundary of the moduli space admits a topological stratification, 
made of loci whose points are curves with a fixed topological type and a prescribed 
assignment of the marked points on each irreducible component. 

The combinatorics of the stable graphs have been investigated in several papers 
in algebraic geometry, for many different purposes (see for instance [ GK) IvOVl) 
lvOV21 IY3| V Our aim with this program is to provide a useful and effective tool 
to generate all the stable graphs of genus g with n unordered marked points up to 
isomorphism, for low values of g and n. 
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We construct an algorithm to generate all the stable graphs of genus g with n 
unordered marked points. Our program uses then the software nauty [M] to elim- 
inate isomorphic graphs from the list of graphs thus created. Since to check that 
two stable graphs are isomorphic is computationally onerous, we try to generate 
a low number of stable graphs, provided that we want at least one for every iso- 
morphism class. The algorithm generates recursively the vectors of genera, number 
of marked points, number of loops, and the adjacency matrix. While it fills these 
data, it checks the stability condition and the condition on the total genus as early 
as possible, in order to minimize the time spent on the branches of the recursion 
that do not lead to stable graphs. Some analysis of the algorithm's performances 
can be seen in Section [H 

Programs for enumerative computations on Aig^n have been implemented in 
both Maple and Macaulay2 f [Fl IY21 lYlj ). Our program can be used, for example, 
to improve the results of |Y3[ Section 5], or to prove combinatorial results on the 
moduli space of pointed stable curves with low genus (cfr. |BMSj . for example 
Corollary 5.3). 

2. Stable modular graphs 

From now on, we fix two natural numbers G and N such that 2G — 2 + N > {). 
For every K E N+, we define I£ = {0, . . . ,K — 1} and E^- to be the symmetric 
group on the set K^. 

Definition 2.1. 

• An undirected multigraph Q is a couple (V, E) with V a finite set of vertices 
and E a finite multiset of edges with elements in y x V/J^2- 

• The multiplicity of the edge (u, w) in E is denoted by mult(ti, w). 

• The total multiplicity oi G, or its number of edges, is \E\: the cardinality of 
E multiset. 

• The degree of a vertex v is degw := 2mult(w,'i;) -I- J^w^^v niult(w, w). 

• A colored undirected multigraph is a multigraph with some additional data 
attached to each vertex. 

Definition 2.2. A stable graph of type {G,N) is a colored undirected multigraph 
Q = {V,E), subject to the following conditions. 

(1) The color of a vertex v is given by a pair of natural numbers (g^, n^). The 
two numbers are called respectively the genus and the number of marked 
points of the vertex v. 

(2) Q is connected. 

(3) Its total genus, defined as J2vev 9^ + l-^l ~ (1^1 ~ equals G. 

(4) Its total number of marked points, defined as X^^jgy '^i" equals N. 

(5) Stability condition: degw -I- n^, > 3 for every vertex v with g^ = 0. 

Notation 2.3. The number degv + is often called the number of half edges 
associated to the vertex v. Condition [s] can be rephrased in: for every vertex v of 
genus 0, its number of half edges is at least 3. 

Two stable graphs G — {V,E,g,n) and Q' = [V ,E' ,g' ,n') are isomorphic if 
there is a bijection /: V ^ V such that: 

• mult(w, w) ~ mult(/(ti), f{w)) for every v,w G V; 

• 9v ^ and n„ = n'^^^^^ for every v gV. 



GENERATING STABLE MODULAR GRAPHS 



3 



Our task is to generate one stable graph for each isomorphism class. 

Remark 2.4. Note that from the definition just given, we are working with an 
unordered set of marked points. The output of the program are the boundary strata 
of the moduli space of stable, genus g curves with n unordered points A^g_„/I]„. 



3. Description of the algorithm 

In this section we describe the general ideas of our algorithm. Let us first intro- 
duce the notation we use in the program. 

Notation 3.1. The set of vertices V will always be i^, so that vertices will be 
identified with natural numbers The multiplicity of the edge between i and 

J will be denoted by a^j : the symmetric matrix a is called the adjacency matrix of 
the stable graph. For convenience, we will denote Ij := aj,j: it is the vector whose 
elements are the number of loops at the vertex j. For simplicity, we will consider 
gj, rij, Ij, ai j to be defined also for i or j outside if, in which case their value is 
always assumed to be 0. 



Remark 3.2. In the following, we assume \V\ > 1 in order not to deal with 
degenerate cases. There are trivially G + 1 stable graphs of type (G, N) with one 
vertex. Indeed, if there is exactly one vertex, the choice of the genus uniquely 



determines the number of loops on it after Definition 2.2 



The program uses recursive functions to generate the data that constitute a stable 
graph. In order, it generates the numbers gj, then the numbers rij, Ij (the diagonal 
part of the matrix a), and finally, row by row, a symmetric matrix representing a. 

When all the data have been generated, it tests that all the conditions of Defi- 
nition |2]2] hold, in particular that the graph is actually connected and satisfies the 
stability conditions. Then it uses the software nauty [M] to check if this graph 
is isomorphic to a previously generated graph. If this is not the case, it adds the 
graph to the list of graphs of genus G with N marked points. 

A priori, for each entry of g, n, and a the program tries to fill that position with 
all the integers. This is of course not possible, indeed it is important to observe 
here that each datum is bounded. From below, a trivial bound is 0, that is, no 
datum can be negative. Instead, a simple upper bound can be given for each entry 
of g by the number G, and for each entry of n by the number TV. For I and a, upper 
bounds are obtained from G using the condition on the total genus (Condition 2.2 1. 

These bounds are coarse: Section [5] will be devoted to proving sharper bounds, 
from above and from below. Also, we will make these bounds dynamical: for 
instance assigning the value go > clearly lowers the bound for gj,j > 0. The 
improvement of these bounds is crucial for the performance of the algorithm. In any 
case, once we know that there are bounds, we are sure that the recursion terminates. 

The algorithm follows this principle: we want to generate the smallest possible 
number of couples of isomorphic stable graphs. To do so, we generalize the idea that 
to generate a vector for every class of vectors of length K modulo permutations, 
the simplest way is to generate vectors whose entries are increasing. The program 
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fills the data row by row in the matrix: 



/ 5o 

no 

lo 



9i 

ni 



9k~i \ 
n-K-i 



(1) 



ai,o 




0-K-l.,K~2 



and generates only matrices whose columns are ordered. Loosely speaking, we mean 
that we are ordering the columns lexicographically, but this requires a bit of care, 
for two reasons: 

• the matrix a needs to be symmetric; in the program we generate only the 
strictly upper triangular part; 

• the diagonal of a need not be considered when deciding if a column is 
greater than or equal to the previous one. 

Therefore, to be precise, we define a relation (order) for adjacent columns. Let 
us call Cj-i and Cj two adjacent columns of the matrix ([l]). They are said to be 
equivalent if cj^i^i — Cj^i for any i ^ {j — 1 + 3, j + 3}. If they are not equivalent, 
denote with jQ the minimum index such that io ^ {j — 1+3, j-|-3} and Cj-i^ip ¥^ Cj,io- 
Then we state the relation Cj-i < Cj if and only if Cj^i^ig < Cj^i^. We do not define 
the relation for non-adjacent columns. We say that the data are ordered when the 
columns are weakly increasing, that is if, for all j, either Cj_i is equivalent to Cj or 

Cj_i <Cj. 

To ensure that the columns are ordered (in the sense we explained before), the 
program keeps track of divisions. We start filling the genus vector g in a non 
decreasing way, and every time a value gj strictly greater than gj-i is assigned, 
we put a division before j. This means that, when assigning the value of Uj, we 
allow the algorithm to start again from instead of nj^i, because the column cj is 
already bigger than the column Cj-i. 

After completing g, we start filling the vector n in such a way that, within two 
divisions, it is non decreasing. Again we introduce a division before j every time 
we assign a value nj strictly greater than nj-i. We follow this procedure also for 
the vector I. 

Finally, we start filling the rows of the matrix a. Here the procedure is a bit 
different. Indeed even if for the purpose of filling the matrix it is enough to deal 
only with the upper triangular part, imposing the conditions that the columns are 
ordered involves also the lower triangular part. A small computation gives that the 
value of j is assigned starting from: 



and we put a division before i if aij > ai-i^j and a division before j if a^j > aij_i. 



max{aij_i, Oi-i,-,} 







if there are divisions before i and j 
if there is a division before i but not before j 
if there is a division before j but not before i 
if there are no divisions before i or j. 
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We cannot conclude immediately that this procedure gives us all possible data 
up to permutations as in the case of a single vector. This is because the transfor- 
mation that the whole matrix undergoes when a permutation is applied is more 
complicated: for the first three rows (the vectors g, n, it just permutes the 
columns, but for the remaining rows, it permutes both rows and columns. Indeed, 
to prove that the procedure of generating only ordered columns does not miss any 
stable graph is the content of the following section. 



We want to prove the following result. 

Proposition 4.1. The algorithm described in the previous section generates at least 
one graph for every isomorphism class of stable graphs. 

From now on, besides G and TV, we also fix the number of vertices K, and focus 
on proving that the algorithm generates at least one graph for every isomorphism 
class of stable graphs with K vertices. 

Notation 4.2. We have decided previously to encode the data of a stable graph 
in a. {K + 3 X K) matrix G := (5, n, /, a) (cfr. ([T])). We denote by A the set of all 
such matrices, and by the set of all {K + 3 x K) matrices that are generated by 
the algorithm described in the previous section. 

We can assume that the graphs generated by the algorithm are stable, since we 
explicitly check connectedness and stability. In other words, we can assume the 
inclusion A4 C A. Hence, in order to prove Proposition |4.1[ we will show that 
every G G ^ is in 7W up to applying a permutation of I£. The idea is to give a 
characterization (Lemma |4.5[ ) of the property of being an element of A4. 

Recall first that the algorithm generates only matrices whose columns are or- 
dered, as described in Section |3] More explicitly, if G = {g,n,l,a) € A, then 
G e M ii and only if: 



Let us call a piece of data gj, rij, Ij, or aij a breaking position if it does not 
satisfy the condition above. Observe that a matrix G £ A has a breaking position 
if and only if G is not an element of A4 . 

We now introduce a total order on the set A of matrices G — {g, n, I, a). If G is 
such a matrix, let v{G) be the vector obtained by juxtaposing the vectors g, n, I 



4. The program generates all graphs 




does not happen, 

^ 9j-i < 9 J , 

gj-i < gj V 7ij-i < Uj , and 

gj-i < gj V < Hj V Ij-i < IjW 



3i' <i:i' ^ {j - A ai/j_i < a.;/ 
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and the rows of the upper triangular part of a. For example, if 

/O 2 0\ 



110 1 




• 111 
1.21 
12.0 
• 



\1 1 

(with the same structure as ([T])), then we define 

z;(G) := (0,0,2,0, 1,1,0,1, 0,0,0,0, 1,1,1, 2,1, 0). 

Definition 4.3. li G,H E A, we write G < H li and only if v{G) is smaller than 
v{H) in the lexicographic order. In this case we say that the matrix G is smaller 
than the matrix H . 

Note that this total order on the set of matrices must not be confused with the 
partial order described in Section [3] From now on we will always refer to the latter 
order on A. 

Remark 4.4. If cr € S^c is a permutation and G = {g,n,l,a) is a graph, then 
we can apply a to the entries of the data of G, obtaining an isomorphic graph. 
The action of ct on G is: {g,n,l,a) — {g' ,n' ,1' , a') where g'^ = gcr{j), = ^iy{j)^ 
I'j = ^a(j) and a[ J = afj[^i)^rj(^jy We denote this new matrix by crG. We write ai_j for 
the element of that corresponds to the transposition of i,j € K_. 



Now we are able to state the characterization we need to prove Proposition 4.1 
Lemma 4.5. Let G £ A; then G €z A4 if and only if G is minimal in the set 

{(T,„i,,G \0<j<K}. 
with respect to the order given in DefinitionU 



Proof. We will prove that G is not minimal if and only if there is a breaking position. 

Assume there is at least one breaking position in G. If there is one in g, n, or 
it is trivial to see that transposing the corresponding index with the previous one 
gives a smaller matrix. If this is not the case, let aij be a breaking position such 
that Oj/j is not a breaking position whenever i' < i (the position {i,j) is the first 
breaking position of its column). We deduce that gj-i = gj, n.j^i = Uj, Ij-i = Ij, 
and that for all i' < i not in {j — l,j}, we have a^'j-i = Oi'j- Let H :— aj-ijG; 
the vectors g, n, and I (the first three rows) coincide in G and H. 

• If J > i, the smallest breaking position is in the upper triangular part of a; 
it is then clear that H < G. 

• If J < i, the smallest breaking position is in the lower triangular part; by 
using the symmetry of the matrix a we again obtain H ^ G (see the right 
part of Figure [I]) . 

Conversely, let j be such that H := aj^ijG -< G. Then consider the first entry 
(reading from left to right) of the vector v{G) that is strictly bigger than v(H). 
This is a breaking position. Notice that if it occurs in the matrix a (equivalently, 
in the last K rows), it is actually the first breaking position of its column. □ 

The proof of Proposition |4.1| follows arguing as in this example. 
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Figure 1. The matrix a when the first breaking position (the 
bullet) is aij with j > i (left) or j < i (right). When transposing 
j — 1 and j, the white and the diagonal-filled entries do not change. 



Example 4.6. Let Gq :— G e Abe the graph of the previous example: 

'0 2 On 



Go 



/ 1 1 1 



• 1 1 1 

1 • 2 1 
1 2 • 



\ 1 2 . 

\ 1 1 • ' 

This graph is stable but not in M because, for example, g2 > gs implies that is 
a breaking position. Thus we apply the permutation <72.3, obtaining the graph 



Gi :— cr2,3Go — 



'0 2^ 
1110 



111 

1 • 1 2 
1 1 • 

• 12 



Now 03,2 is a breaking position; applying ai^2, we obtain 



G2 :— cri,2Gi — 



/ 1 1 1 

0000 

• 1 1 1 
1 • 1 
1 1 • 2 



\ 1 1 • 2 

V 1 2 • ' 

This introduces a new breaking position at 03^1, so we apply the transposition cfoX- 



2n 
1110 



• 1 1 

1 • 1 1 



-< G-2 



G3 :— (To, 102 — 

\ i 1 i 2 
\ 1 2 • ' 

The graph G3 is finally in Ai and indeed no transposition can make it smaller. 



Proof of Proposition \4-l\ Recall that we have to prove that for every G £ A, there 
is a permutation a e such that aG € M. 

So, let Go = G e yl. If G e A^, then we are done; otherwise, G does not 
satisfy the condition of Lemma 4.5 hence there is a transposition crj_i.j such that 
Gi — aj-ijGo -< Gq. 

The iteration of this process comes to an end (that is, we arrive to a matrix in 
Ai) since the set 

{aG\ae^K} 

is finite. □ 



8 



STEFANO MAGGIOLO AND NICOLA PAGANI 



5. Description of the ranges 

In Section [3] we have introduced the algorithm, by describing the divisions. In 
this section we introduce accurate ranges for the possible values of g, n, I and a. 



We will deduce from the conditions of Definition 2.2 some other necessary condi- 
tions that can be checked before the graph is defined in its entirety. More precisely, 
every single datum is assigned trying all the possibilities within a range that de- 
pends upon the values of G and N, and upon the values of the data that have 
already been filled. The conditions we describe in the following are not the only 
ones possible; we tried other possibilities, but heuristically the others we tried did 
not give any improvement. 

The order in which we assign the value of the data is g, n, I, and finally the 
upper triangular part of a row after row. 

Notation 5.1. Suppose we are assigning the i-th value of one of the vectors g, n 
or I, or the (i,j)-th value of a. We define the following derived variables e™'^^, c 
and pi that depend upon the values that have already been assigned to g, n, I, a. 

We let e™^^ be the maximum number of edges that could be introduced in the 
subsequent iterations of the recursion, and c be the number of couples of (different) 
vertices already connected by an edge. We let pi be the number of vertices z 
to which the algorithm has assigned gz — 0. Note that the final value of pi is 
determined when the first genus greater than is assigned, in particular the final 
value of pi is determined at the end of the assignment of the values to the vector 
g. On the other hand, c starts to change its value only when the matrix a begins 
to be filled. 

After the assignment of the i-th value, the derived values e™'*'^, c and pi are then 
updated according to the assignment itself. 

(2) 

Notation 5.2. When deciding g, n, or /, we let be the minimum between 2 
and the number of half edges already assigned to the i-th vertex. This is justified 
by the fact that we know that, when we will fill the matrix a, we will increase by 
one the number of half edges at the vertex i in order to connect it to the rest of the 

(2) 

graph. Hence, whenever gi = 0, is the number of stabilizing half edges at the 
vertex i: one half edge is needed to connect the vertex to the rest of the graph, and 
then at least two more half edges are needed to stabilize the vertex. When deciding 
ttij, it is also useful to have defined hi, the total number of half edges that hit the 
i-th vertex. Finally, we define 









E ' 




i' <i 






7V(2) 






. (2) 








i' < i 
















:= E"*'J' 




i' <i 




i'<iVj'<j 



We are now ready to describe the ranges in which the data can vary. We study 
subsequently the cases of g, n, I and a, thus following the order of the recursions 
of our algorithm. Each range is described by presenting a first list of general 
constraints on the parameters and then by presenting a second list containing the 
actual ranges in the last line. 
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5.1. Range for gi. When the algorithm is deciding the value of gi, we have the 
following situation: 

• e""'"'' ^G~G, + K -Ihy Condition H 

• amongst the e™''^ edges, there are necessarily K — 1 non-loop edges (to 
connect the graph) ; these K — 1 edges give one half edge for each vertex, 
whereas we can choose arbitrarily where to send the other K — 2 half edges; 
conversely, the 2(e™'^'^ — K + 1) half edges of the remaining edges can be 
associated to any vertex; therefore, the maximum number of half edges (not 
counting those that are needed to connect the graph) is 2e™'*'' — K + N = 
2{G-Gi)+K -2 + N- 

• we need 2pi half edges to stabilize the genus vertices, since one half edge 
comes for free from the connection of the graph. 

We use the following conditions to limit the choices we have for gi: 

(1) since g is the first vector to be generated, there is no division before i, hence 

9i > 9i^u 

remember that gj = whenever j ^ K;, 

(2) we need at least K— \ non-loop edges, hence (using the fact that X]j>i 9j — 
{K-i)gi) 

^ G - G, - {K ^ i)g, +K -1>K ~1 
^{K- i)g, <G-G,: 

(3) in order to stabilize the pi vertices of genus (using the fact that one 
stabilizing half edge comes for free by connection) we must have 

2pi < 2e"^'' - K + N 

^2pi<G-G,-{K- i)g, -K + N 
^{K~ i)gi <G-G,- K + N -2pi. 

5.2. Range for rii. When deciding n^, we have the following situation: 

• as before, e™'^^ — G — Gk + K — 1>K — 1, and the maximum number of 
half edges stiU to be assigned is 26"^^=^ - K + N - N, ~ n, ^ 2{G - Gk) + 
K -2 + N - Ni-n,; 

(2) (2) 

• we need 2pi — — nl half edges to stabilize the first pi vertices; 

(2) (2) 

• if (7i = 0, we need 2{i + 1) — more half edges to stabilize the first 
i + 1 vertices. 

The following conditions define then the ranges for the possible choices for n^: 

(1) if there is not a division before i (that is, if gi — gi-i), then we require 
rii > rii-i] otherwise, just rii > 0; 

(2) we cannot assign more than TV marked points, hence (where we treat the 
case of = in a special way) 

Ni + Ui < N 

^ Ui < N - Ni 

(pi — i)ni < N — Ni li moreover gi = 0. 
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(3) if (7i = 0, for the purpose of stabilizing the first i + 1 curves we cannot use 
marked points anymore, therefore we have 

2{i + 1) - A^f ^ - nf ^ < (2(G - Gk) + K ~ 2) 

^nf '> =min(2,n,) > -{2{G - Gk) + K - 2) + {2{i + 1) - N^^'') 

{impossible if RHS > 2 
Ui > RHS otherwise. 

5.3. Range for li. When deciding li, this is the situation: 

• e™'^'^ = G — Gk — Li — li + K— 1>K~1, and the maximum number of 
half edges still to assign is 2e™^'^ - A' = 2(G - Gk - Li - k) + K - 2; 

The conditions on li are then the following: 

(1) if there is not a division before i, then we require li > U-i; otherwise, just 
k > 0; 

(2) we need at least K — 1 non-loop edges, hence 

=>G-GK~L,-k+K-l>K-l 
^ h < G — Gk — Li ; 

(3) let z be the index of the genus vertex with the least number of stabilizing 
half edges such that z < i; it already has + 2lz half edges, but we cannot 
use loops anymore to stabilize it; hence, 

max(0, 2-71^ - 2h) < G - Gk - Li - k + K - I 

^ k < G - Gk - L, + K - 3 + + 21^ 

(4) assume gi = 0; if li > 0, we are adding to the i-th vertex 2 — nl' stabilizing 
half edges, and to stabilize the pi genus vertices, we need to have 

2pi - iV(2) _ (2 - nf ^) < 26°"^'^ - K 

^ 2pi - iV(2) - (2 - nf^) max(0, 2 ~ rrii) <2{G - Gk - L, - k + K - I) - K 

^ 2k < 2{G - Gk - L,) +K + iV^^) _ „(2) _ 2p, . 

(5) assume gi = 0; after deciding k, we still have e"^^^ edges to place, and 
each of them can contribute with one half edge to the stabilization of the 
i-th vertex; moreover, one of these half edges is already counted for the 
stabilization; hence 

n, + 2k + (e"'''^ - 1) > 2 

^ + 2k + G - Gk - Li - k + K - I - 1 > 2 
^ li>A- n.i- G + Gk + L^- K . 

5.4. Range for ai_j. When deciding a.^j, this is the situation: 

• earlier in Notation |5.2[ we observed that for the purpose of filling the vectors 
g,n and / we could consider a genus vertex stabilized when it had at least 
two half-edges (since the graph is going to be connected eventually) . When 
assigning the values of a, the stability condition goes back to its original 
meaning, i.e. each vertex has at least 3 half edges. 

• e'""'^ = G-Gk-Lk- A,,j +K-1; 

• we have already placed edges between c couples of different vertices; 
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Here are the constraints that aij must satisfy: 

(1) if there is not a division before i, then we require Uij > flj-ij-; otherwise, 
just ttij > 0; 

(2) if there is not a division before j, then we require aij > ajj-i; 

(3) we need at least K—2 — c {if positive) edges to connect the graph, because if 
ai_j > 0, c will increase by 1 (this estimate could be very poor, but enforcing 
the connectedness condition in its entirety before completing the graph is 
too slow), hence: 

e™'*'^ - ttij > max(0, K-2-c) 

aij < G - Gk - Lk - Aij +K -1- max(0, K -2-c); 

(4) Uij contributes with at most max(0, 3 — /ij) H-max(0, 3 — /ij) stabilizing half 
edges; hence, to stabilize the pi genus vertices, we need 

3pi - ^ min(3, rij + 2li) - (max(0, 3 - hi) + max(0, 3 - hj)) < 2{e"'^'' - a^j) 

3pi — ^ min(3, rij + 2li) — (max(0, 3 — ftj) + max(0, 3 — hj)) < 

< 2(G ~Gk-Lk- A,,j +K -1- a,;j ) 
=^ 2o,,j < 2{G -Gk- Lk- A,,^ + K - 1) - 3pi + 

+ min(3, rij + 21^) + max(0, 3 — /ij) + max(0, 3 — hj) . 

(5) if j = K — 1 (that is, if this is the last chance to add half edges to the i-th 

vertex) , then wc add enough edges from z to ii" — 1 in order to stabilize the 
vertex i; moreover, if up to now we did not place any non-loop edge on the 
vertex i, we impose ai^K-i > 0. 

0'i,K-i > if = for all 1 < J < ii" — 1, 

ai^K-i >3-hi if fifj = 0. 

6. Performance 

The complexity of the problem we are trying to solve is intrinsically higher 

than polynomial, because already the amount of data to generate increases (at 
least) exponentially with the genera and the number of marked points. We also 
observed an exponential growth of the ratio between the time required to solve 
an instance of the problem and the number of graphs generated. Anyway, our 
program is specifically designed to attack the problem of stable graphs, and it can 
be expected to perform better than any general method to generate graphs applied 
to our situation. 

We present here some of the results obtained when testing our program on an 
Intel® Core''''^2 Quad Processor Q9450 at 2.66 GHz. The version we tested is not 
designed for parallel processing, hence it used only one of the four cores available. 

However, when computing a specific graph, the program needs to keep in the 
memory only the graphs with the same values in the vectors g, n, I: memory 
usage becomes therefore negligible. Moreover this shows that we can assign the 
computations of stable graphs with prescribed <?, n, I to different cores or CPUs, 
thus having a highly parallelized implementation of the program. 
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Figure 2. Time needed to compute all stable graphs of type (G, N). 



G 


N 


Time (s) 


# stable graphs 


Duplicates (%) 





18 


392 


847,511 


54.9 


1 


14 


539 


1,832,119 


41.3 


2 


10 


147 


1,282,008 


30.9 


3 


7 


117 


1,280,752 


29.9 


4 


5 


459 


2,543,211 


40.1 


5 


3 


606 


2,575,193 


54.7 


6 


1 


226 


962,172 


70.6 


7 





681 


1,281,678 


85.6 



Table 1. For small G, the maximum N such that all stable graphs 
of type (G, N) can be computed in less than 15 minutes. In the 
last column we show the ratio of duplicated graphs among the total 
number of those created by our generation algorithm. 



In Table [T] we list, for each genus G, the maximum number of marked points N 
for which we can compute all the stable graphs of type (G, N) under 15 minutes. 

In Figure |2] we show all the couples (G,iV) that we computed against the time 
needed; the lines connect the results referring to the same genus. From this plot it 
seems that, for fixed G, the required time increases exponentially with N. However, 
we believe that in the long run the behaviour will be worse than exponential. This is 
suggested also by the fact that the ratio of non-isomorphic stable graphs over those 
created by our generation algorithm tends to zero as G and N grow (see Figure |3]). 

More benchmarks and up-to-date computed results are available at boundary's 
webpage, |http : / /people . sissa. it/~maggiolo/boundary/[ 
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Figure 3. Ratio of duplicated graphs among the total number of 
those created by our generation algorithm. 
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